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One-loop amplitudes of gluons in M = 4 gauge theory can be written as linear com- 
binations of known scalar box integrals with coefficients that are rational functions. In 
this paper we show how to use generalized unitarity to basically read off the coefficients. 
The generalized unitarity cuts we use are quadruple cuts. These can be directly applied to 
the computation of four-mass scalar integral coefficients, and we explicitly present results 
in next-to-next-to-MHV amplitudes. For scalar box functions with at least one massless 

external leg we show that by doing the computation in signature ( h+) the coefficients 

can also be obtained from quadruple cuts, which are not useful in Minkowski signature. 
As examples, we reproduce the coefficients of some one-, two-, and three-mass scalar box 
integrals of the seven-gluon next-to-MHV amplitude, and we compute several classes of 
three-mass and two-mass-hard coefficients of next-to-MHV amplitudes to all multiplicities. 
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1. Introduction 



One-loop amplitudes of gluons in supersymmetric field theories are four-dimensional 
cut constructible |l],0. That is to say, they are completely determined by their unitarity 
cuts. 

In the particular case of Af = 4 gauge theories, the amplitudes can be written as a sum 
over scalar box integrals with rational functions as coefficients ]TJ||]. Scalar box integrals 
are one-loop Feynman integrals in a scalar field theory with four external legs and four 
propagators. The scalar in the loop is massless, and the ones in the external legs could be 
massive. This gives four families classified by the number of massive legs. These are called 
one-, two-, three-, and four-mass scalar box integrals. Given that all scalar box integrals 
are known explicitly ||, the task of computing a given one-loop amplitude is reduced to 
that of finding the coefficients. 

These coefficients were computed for all maximally helicity violating amplitudes in 
Jlj], for next-to-MHV six-gluon amplitudes in ||, for next-to-MHV seven-gluon amplitude 
with like helicity gluons adjacent in 0] and for all helicity configurations in ||. 

One common feature of all these computations is that they are based on the use of 
unitarity cuts. The basic idea is to compute the discontinuity of the amplitude across a 
given branch cut in the kinematical space of invariants by adding all Feynman graphs with 
the same cut. The discontinuity is obtained by "cutting" two propagators. This sum of 
cut Feynman graphs reduces to the product of two tree-level amplitudes integrated over 
the Lorentz-invariant phase space of cut propagators. On the other hand, each of the 
scalar box integrals are also Feynman graphs whose discontinuity can also be computed 
by cutting two propagators. Therefore, combining the two different ways of computing the 
same discontinuity, one can get information about the coefficients. 

One complication of all these approaches is that there are several scalar box integrals 
sharing a given branch cut. Therefore, one has to disentangle the information of several 
coefficients at once. This is done in [0],[2|,[)| by using reduction techniques an d in ||[|] 

by using the holomorphic anomaly ||10|], which affects the action of differential operators 



that test localizations in twistor space 11 



In this paper we present a different way of computing the coefficients, using generalized 
unitarity. Even though several scalar box integrals can share the same branch cut, it 
turns out that their leading singularity is unique. For a detailed treatment of generalized 



unitarity and leading singularities of Feynman graphs, see Sections 2.9 and 2.2 of ||12||. 



Therefore, by studying the discontinuity associated to the leading singularity, which we 
denote by A^g, we can basically read off a given coefficient. 

For a given Feynman graph, A^s is obtained by cutting all propagators. In |T3[ the 
leading singularity of a three-mass triangle graph, whose discontinuity is computed by a 
triple cut, was used to compute its contribution to the e + e~ — > ( i y*,Z) — > qqgg one-loop 
amplitude. In this paper, we use a quadruple cut to compute coefficients of four-mass box 
integrals in one-loop M = 4 gluon amplitudes. We then turn our attention to scalar box 
integrals with at least one massless leg. In Minkowski space, Als for these box integrals, 
which is still a quadruple cut, does not give information about the coefficients, as we 
explain in section 2. One has to use at most triple cuts.0 One disadvantage of triple cuts is 
that there might be several box integrals with the same singularity, and therefore several 
unknown coefficients will show up at once. 

We find a way out by going to signature ( h+).@ We show that in this signature 

all scalar box integral coefficients can be computed by a quadruple cut. This allows us to 
also read off the coefficients directly, since again only one coefficient contributes in a given 
cut. The quadruple cut integral is completely localized by the four delta functions of the 
cut propagators. It reduces to a product of four tree-level amplitudes, which can be easily 
computed by MHV diagrams [p~^| , pT5|JT^ ] . This is a very simple way of computing any given 
coefficient in any one-loop M = 4 gluon amplitude. 

We illustrate this procedure by computing several coefficients of one-, two-, three- and 
four-mass box scalar integrals. For the first three cases we give examples for next-to-MHV 
amplitudes, including some results to all multiplicities, and for the last we give examples 
for an eight-gluon next-to-next-to-MHV (NNMHV) amplitude. 

Motivated by the introduction of twistor string theory ||11|| , there has recently been 



interest in the twistor-space localization of gauge theory amplitudes [jT4] , |T71jn| ■ In partic- 
ular, for supersymmetric gauge theories, the coefficients of box, triangle and bubble scalar 
integrals also exhibit simple twistor space structure [f^ilf. At one- loop, this structure can 
be understood from the MHV diagram formulation ||f 9||. 



1 Triple cuts were used in @ to conclude that some coefficients in one-loop M = 4 gluon 
amplitudes must vanish for special choices of helicity. 

2 Alternatively, it is perhaps more natural to drop the restriction of real momenta and work 
with the complexified Lorentz group. We thank E. Witten for pointing this out. 
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Since no four-mass coefficient has previously been presented in the literature, we 
consider its twistor space support. This is done by studying various differential operators 
acting on unitarity cuts. 

This paper is organized as follows: In section 2, we discuss the generalized unitarity 
method in general and discuss how one can use quadruple cuts for all scalar box integrals in 

signature ( (-+)• In section 3, we present examples of each type of scalar box integral, 

starting with the four-mass, and then the three-, two-, and one-mass. In section 4, we 
give several examples of infinite classes of coefficients in next-to-MHV amplitudes. In the 
appendices, we consider in detail the discontinuity of the four-mass scalar box integral 
associated to branch cuts in all possible channels and use this to get consistency checks on 
the new coefficient obtained in section 3. 



1.1. Preliminaries 

One-loop Af = 4 amplitudes of gluons can be written as a linear combination of 
scalar box integrals with rational coefficients In this paper, we concentrate on 

the leading-color contribution, which is the part of the full amplitude proportional to 
iVTr (T ai . . .T a ™). 

We write this schematically as follows: 

A n:1 = J2 (bl lm + cl 2m e + dl 2m h + gl 3m + // 4m ) . (1.1) 

The integrals are defined in dimensional regularization as 

2- f d 4 ~ 2e £ 1 
I 4 = -i(4n) e y ( 27 r) 4 - 2e i\i - Ktf{l -K x - K 2 ) 2 (£ + K 4 ) 2 ' ^ 

The external momenta Ki are taken to be outgoing and are given by the sum of the mo- 
menta of consecutive external gluons, as shown in Figure 1. The labels lm, 2m, 3m, 4m 
refer to the number of legs K n such that K% ^ 0, or equivalently, the number of vertices 
in the box with more than one external gluon. For two-mass box integrals, there are two 
inequivalent arrangements of massive legs. Either they are adjacent (J 2m or they are di- 
agonally opposite (J 2m e ). All these integrals are UV finite but suffer from IR divergences, 
except for 7 4m which is finite. 
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Fig. 1: Scalar box integrals, (a) The outgoing external momenta at each of the 
vertices are K\ , K2 , K3 , K4 , defined to correspond to sums of the momenta of gluons 
in the exact orientation shown, (b) One-mass I^a- (c) Two-mass "easy" 1%™-,?- 
(d) Two-mass "hard" 7, 
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From here on, we will drop the dimensional regularization parameter e, because we 
will only deal with cuts that are finite. Moreover, we work in four-dimensional Minkowski 
space. 

In the literature, it is common to write the amplitude in terms of scalar box functions. 
These functions are given in terms of the scalar box integrals as follows. 



rplm Tp2m e rp2m h 

jim _ 2 n:i j 2m e = _2 n ' r ' 1 I 2rn h = —2 n:r;i 



;i.3) 



4:i ~ K 2 K 2 ' 4:r ' i K 2 K 2 — K 2 K 2 ' 4:r;i ~~ K 2 K 2 ' 

JX A1 JX l2 iY 41 iY 12 ft l A 3 i M2 i Ml 

rp3m rpAm 

jSm _ _r> n:r:r';i j-4 m _ Q n:r:r' :r" ;i 

A 4::r:r';i — z K 2 K 2 K 2 K 2 ' 4:r:r' :r" ;i ~ K 2 K 2 _ ' 

- n -41- n -12 ^1^3 -"-41-^12^ 

Here, and throughout the paper, we define 

K = K + K 

p= V / 1-2A 1 -2A 2 + (A 1 -A 2 ) 2 , 

A _ zi^L (i.4) 

At 



/^ 2 



Then we can alternatively write (|1 . 1| ) as a linear combination of scalar box functions 



A n:1 = i bFlm + ci?2m 6 + dF2m h + 9 F3m + fF 4 ™) ■ ( L5 ) 

Each way of writing the amplitude has its own advantages and disadvantages. In ( |1 . 1| ) , 
all coefficients are rational, but their twistor space support is not simple. In ( |1.5| ), the 
coefficient of the four-mass box function is not rational, for it contains a square root, but 
all coefficients have simple twistor space structure. For a discussion of the localization 
in twistor space of the four-mass scalar box function coefficient, see the appendices. For 
the rest of the body of the paper, we will work mainly with scalar box integrals and their 
coefficients as formulated in (|1 . 1|) . 



2. Generalized Unitarity and Quadruple Cuts 

One-loop amplitudes in field theory have several singularities as complex functions of 
the kinematical invariants. In M = 4 gauge theory, the singularities can only be those 
of the scalar box integrals (|1 . 2| ) and of the coefficients in ( |1 . 1|) . Since the coefficients are 
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rational functions, they are not affected by branch cut singularities. Therefore one can get 
information about them by studying the discontinuities of the amplitude across the cuts. 

In fact, most of the techniques for computing the coefficients efficiently are based on 
studying unitarity cuts . The basic idea is to consider the branch cut singularity 

of the amplitude in a given channel. The discontinuity across this branch cut can be 
computed by a cut integral on both sides of the equation 

A n:1 = + cl 2m e + dl 2m h + gl 3m + fl 4m ^j . (2.1) 

On the left hand side, one cuts two propagators of all Feynman integrals participating in 
this channel, while on the right hand side one cuts two propagators of scalar box functions. 

By "cutting propagators" we mean the following. In Minkowski space, propagators in 
one loop integrals are defined by using Feynman's ie prescription, i.e., l/(-P 2 + ie). This is 
equal to the principal value of 1/P 2 plus i5^ + \P 2 ), where (+) indicates a restriction to the 
future light-cone. Cutting a propagator means removing the principal part, i.e., replacing 
the propagator by i5^ + \P 2 ). 

The sum over cut Feynman integrals on the left of (|2.1| ) becomes an integral over a 
Lorentz invariant phase space of the product of two tree-level amplitudes. To be more 
explicit, consider the cut in the (i, i + ...^-channel, 

C = J d\i A tree (£i, i, £ 2 ) A tree (-£ 2 ,j + 1, - 1, —£\), (2.2) 

where d\i is the Lorentz invariant phase space measure for (£i, £ 2 ). It is given explicitly by 

dix = ^\l\) S^ih+h-P^), (2.3) 

where Pjj is the sum of the momenta of gluons from i to j. 

If we denote the discontinuity across the branch cut of a given scalar box integral I 
by AI, then we have the following equation. 

f dfxA^i^i, j, e 2 )A t ^(-£ 2 J + 1 -h) = 

J (2.4) 

(bAI lm + cAI 2m e + dAI 2m h + gAI 3m + fAI 4m ^ . 

As discussed in the introduction, these equations have been used to get the coefficients for 
MHV amplitudes and six- and seven- gluon next-to-MHV amplitudes 0,0,^^,0. 
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However, extracting the coefficients is not a simple task in general. The main reason is 
that several scalar box integrals share a given cut, and therefore their unknown coefficients 
enter in the equation at the same time. 

Several approaches exist to extract the coefficients. One is based on reduction tech- 
niques that allow writing the integrand of ( [2.2| ) as a sum of terms that have the structure 
of cut scalar box functions [0,0. Another method |||[] uses operators that test localization 
in twistor space to get rational functions on both sides of the discontinuity of ( |2~7i] ) and 
compare the pole structure. 

It is the aim of this section to use higher order singularities to reduce the number of 
scalar box functions that enter in the generalization of ( |2.4|) . Ideally, we would like to find 
only one scalar box integral on the right hand side of ( |2.4| ). 

This can easily be done for the four-mass scalar box integral. Consider the disconti- 
nuity associated to its leading singularity, Als- As mentioned in the introduction, this is 
computed from the integral 



r4m 



dH {l 2 + ie)((£ - i^i) 2 + ie)((£ - K x - K 2 ) 2 + ie)((£ + K 4 ) 2 + ie) (2 ' 5) 



by cutting all four propagators: 
A LS I 4m = j d^£5 {+) {£ 2 ) c^ + )((£-Ki) 2 ) 5 {+ \{£ - K x - K 2 ) 2 ) S i+) ((£ + K 4 ) 2 ). (2.6) 




m+1 



Fig. 2: A quadruple cut diagram. Momenta in the cut propagators flows clockwise 
and external momenta are taken outgoing. The tree-level amplitude A* ree , for 
example, has external momenta i + 1, j, £2,£i- 
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a tree a tree a tree a tree ?a r4m 
A l A 2 A 3 A 4 = J^LSJ- , 

where A^ ree is the tree-level amplitude at the corner with total external momentum K,, 
See fig. 2. 



Note that since there are four delta functions, and £ is a vector in four dimensions, 
then the integral is localized and equals a jacobian J = {A.K\ x K\ 2 p)~ x . Moreover, the 



where the sum is over the possible spins J of internal particles and the solution set S of ts 
for the localization equations, |«S| is the number of these solutions, and rij is the number 
of particles of spin J. 

This equation gives the coefficient / of any four-mass box function in any one-loop 
amplitude, as we discuss in more detail in section 2.1. 



same jacobian appears on both sides of (|2.7| ), and we find that 




(2.8) 
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The natural question at this point is whether the same thing can be done for other 
scalar box integral coefficients. It turns out that in Minkowski space this is not possible. 
Recall that all other scalar box integrals contain at least one massless vertex. This means 
that at least one of the four tree-level amplitudes in the quadruple cut must be a three- 
gluon amplitude. One problem is that all three-gluon amplitudes vanish on-shell, as will 
become clear in the discussion below. 

For definiteness, consider the three-mass scalar box integral in Figure 3. Let us denote 
the momenta of the three-gluon vertex by £±, £4 and p. We have to impose that all three 
vectors be lightlike and that £4 = £\ +p. From this it is easy to see that necessary conditions 
for all three vectors to be lightlike are £\ ■ p = £4 ■ p = £\ ■ £4 = 0. In Minkowski signature 
one can go to a frame where p = (E, 0, 0, E). Then it is easy to see from l\ • p = that 
£\ = ap. Likewise, £4 = Pp. This means that all three gluons are collinear. Momentum 
conservation then implies that f3 = a + 1. 

On the other hand, from £2 = £\ — K\ and £\ = we can compute a, and from 
£3 = ^4 + K3 we get (3. They are given by 

K 2 K 2 
2K\ ■ p 2K 3 ■ p 

Now we can ask if the equation f3 = a + 1 is satisfied. It is not difficult to see that this 
equation is equivalent to K 2 K 2 — {K\ + K^f'iKx +p) 2 = 0. Note that the left hand side is 
precisely the denominator in (|1.3| ). This means that this singularity is a pole, as opposed 
to the usual branch cut singularity. The discontinuity associated to it is a delta function,! 
reflecting the fact that only three of the four delta functions suffice to perform the integral. 
This implies that the quadruple cut equation ( p.7|) is in this case given by 



J dfi A\ ree A l 2 ee A^ ee A l l ee = g 6(KfK$ - {K x + K 2 f(K x + p) 2 ). (2.10) 

Note that the left hand side is trivially zero, for A4 is a three-gluon tree-level amplitude 
on-shell which vanishes in Minkowski space. Then we find that 

g 8{K 2 x Kl - (Kt + K 2 f{K x + p) 2 ) = 0, (2.11) 



3 In one real variable, the statement that the discontinuity associated to a pole is a delta 
function is expressed in the familiar relation — ; }-r- = 2m8(x). 
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and therefore the coefficient must have a zero at the support of the delta function. In- 
deed, recall that the relation between the box integral coefficient 'g and the box function 
coefficient g is given by g = ^(KfR^ - (K x + K 2 ) 2 (K 1 + p) 2 )g. 

This proves that the information about g disappears from the equation, and the 
quadruple cut cannot be used to learn anything about g, except for the presence of the 
factor \{K\Kl - {K x + K 2 f{K x + pf). 

If we insist on using Minkowski signature, we have to study triple cuts.ll The disad- 
vantage is that in most cases, more than one box integral will share the same singularity. 
The situation is better than with the usual cut in ( |2.4| ), but it is still complicated in most 
cases. 

2.1. Quadruple Cut in Signature ( h+) 

The main source of the problem for scalar box integrals with at least one massless leg 
is the presence of the three-gluon vertex. The fact that this always vanishes in Minkowski 
space made impossible the extraction of their coefficients from the quadruple cut. 

As explained carefully in section 3 of [[TTj, three-gluon amplitudes with helicities (+ + 
— ) and ( h) do not necessarily vanish in other signatures. In particular, motivated by 



the study of a string theory with target twistor space [11], Witten considered the three- 
gluon amplitude in signature ( h+). It turns out that this is precisely what is needed 

in order to use the quadruple cuts to read off all coefficients. 



We now review the analysis done in ||11||. In four dimensions we can write a null 



vector as a bispinor, p a a = X a Xa- The inner product of two null vectors p and q is given by 
2p ■ q = (X p Xq)[X p X q ], where the brackets represent the natural inner products of spinors 
of positive and negative chirality. 

A tree- level three-gluon amplitude with helicity (++— ) or ( h) is given respectively 

by pr 



A3 (x r x p][ x q A r f A " {p ' r } " (K x p) (x q x r y (2 ' 12) 

In Minkowski space and for real momenta, X p and X p are complex but not independent, 
X p = ±A p . Therefore, when p ■ q = 0, this means that both (p q) and [p q] vanish. This 



implies that both amplitudes in (|2.12|) also vanish. 



4 To be precise, all unitarity cut analyses require complexified momenta. We thank L. Dixon 
for this clarification. 
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On the other hand, in signature ( h+) and for real momenta, both A p and X p are 

real and independent. Therefore, 2p ■ q = (X p X q )[X p X q ] = has two solutions. Either 
(X p X q ) = or [A p A,] = 0. 

Since momentum conservation implies that p-q = p- r = q- r = 0, we find that if 
(A p X q ) = 0, then (X p X r ) =0 and (A r X q ) = must also hold. Therefore all three A's are 
proportional. Likewise if we choose [A p X q ] = 0, then all three A's are proportional. 

Now it is clear that if we are faced with a (+ H — ) tree amplitude, we should choose 
all A's to be proportional, and then the corresponding amplitude in ( |2.12| ) will not vanish. 
Likewise, if we are faced with a ( h), we should choose all A's to be proportional. 

Having solved the problem of the vanishing of the amplitudes, we have to deal with 

the meaning of a "unitarity cut" in signature ( We certainly cannot offer a formal 

theory here, but some interesting developments have appeared in the literature about field 
theories in ( h+) signature. For a recent review, see [21]. 

Here we take a more operational approach in order to compute the coefficients. Since 
the coefficients are written in terms of invariant products of spinors, once they are com- 
puted they can be used in Minkowski space. 

Consider the quadruple cut measure in Minkowski space. Recall that each delta 
function is actually given by 

£(+)(P 2 ) = $(E P )5(P 2 ), (2.13) 

where Ep is the zeroth component of P M . Both sides of ( |2.7|) contain the same measure 
and in particular the same factors of $(x). Since the integral is localized, one can drop the 
$(x) factors on both sides. We then take as our definition of a quadruple cut in signature 
( h+) the following: 

/ dH 5(£ 2 ) 8{(t - K x f) 8{{i -K 1 - K 2 ) 2 ) 5((£ + K^A^A^A^A^ = 
J , (2-14) 
g J dH 5(£ 2 ) 5((£ - K x f) 5((£ -K x - K 2 ) 2 ) S((£ + K,) 2 ). 

Here we use ^ to represent a generic coefficient in ( |2.1|) with at least one massless leg. 

In ( |2.14p the delta functions are ordinary Dirac delta functions. One can solve for 
£ a a, taking into account that solutions for which one or more of the tree- level amplitudes 
vanish do not contribute. Therefore we find that the coefficient must be equal to 

g = j^J2 ^j^i ree ^2 ree ^3 ree ^4 ee , (2-15) 

' ' S,J 

where again the sum is over all solutions. 

In the rest of the paper, we illustrate this procedure in detail. 
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3. Examples 

The previous section demonstrated that coefficients of box integrals may be computed 
from the formula 

S - = |4[ E ^J< ee ^2 ree ^3 ree < ee , (3-1) 

where J is the spin of a particle in the M = 4 multiplet and S is the set of all solutions of 
the on-shell conditions for the internal lines, 

5 = {£\£ 2 = 0, (£-K 1 ) 2 = 0, {i-K x - K 2 f = 0, (£ + K 4 ) 2 = 0}. (3.2) 

In this section, we present a variety of applications of the formula ( |3.1| ). 
The explicit covariant solution of ( |3.2|) for the vector £ is given by: 

£ = PxKt + (3 2 K 2 + f3 3 K 4 + (3 4 P- 
P — e K v K p K a 

6 = ^2 K^2 ■ ■ ^4)(2(^i ■ K 2 ) + Kl) + K 2 (K 2 ■ K 4 )) 
+(2(K 1 ■ K 2 f - (K 2 + {K x ■ K 4 ))K 2 + {K x ■ K 2 ){K 2 + [K 2 ■ K 4 )))K 2 ] , 

& = ^2 K^i • KiW&i ■ K 2)(Ki ■ K 4 ) + (K, ■ K 4 )K 2 - K 2 (K 2 ■ K 4 )) (3 . 3) 

-(-(ifi ■ K 2 )(K! ■ K 4 ) + Kl({K x ■ K 2 ) + K 2 + (K 2 ■ K 4 )))K 2 } , 
& = ^2 K^i • K 2){~{K X ■ K 4 )K 2 + K 2 (K 2 ■ K 4 )) 

-{Kx ■ K 2 ) 2 (2(K 1 ■ K 4 ) + Kl) + K 2 K 2 {{K X ■ K 4 ) + (K 2 ■ K 4 ) + K 2 )] , 
±K 2 2 K 2 lP 



I3 4 = 



AP 2 



In particular, |«S| = 2, and the two solutions S+, S- are related by a change of sign of p, 
which appears in (3 4 . Since there is a sum over the two solutions S+ and S-, the result for 
a a is seen to be rational, as it must be @.@ 

We begin with the four-mass box integral, where the solution ( |3.3| ) can be used directly 
in the formula ( |3.1| ) .i 



5 For one-, two-, and three-mass coefficients, p turns out to be rational, so each solution is 
individually rational. 

6 The formula ( |3.3D is implicitly written in signature ( 1~+)- To perform a calculation with 

external momenta in Minkowski space, one must Wick-rotate e M „ pCT — > itfivpa- 
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For all other box integrals, we must work in the signature ( 1~+)> as described 

in section 2. As we have seen, the helicity configuration determines how to solve for the 
spinor components of the cut propagators. 

It is worth noting that for all these cases, there are still just the two solutions ( |3.3[ ); in 
fact, a given solution determines whether the holomorphic or antiholomorphic spinor s are 
proportional at each three-gluon vertex. This relation is illustrated in Figure 4, for each 
family of box functions with a three-gluon vertex. 
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Fig. 4: Possible helicity assignments at three-gluon vertices, as derived from the 
two solutions S+ , S- given in ( |3.2| ) . The assignment + H — means that the solution 
dictates that the holomorphic spinors A are all proportional at this vertex; the 

assignment h means that the solution dictates that the antiholomorphic spinors 

A are proportional. The three helicities can then be distributed freely on the three 
legs at that vertex. 



In practice, however, it is useful to solve for spinor components directly, rather than 
using ( |3.3|) as written. This is because the four tree-level amplitudes in ( |3.1|) can be easily 



computed from Parke- Taylor formulas []20] or MHV diagrams [fP4"| , |T5|| , which involve spinor 
products. 

In the three-mass example, we demonstrate the judicious solution of spinor compo- 
nents. In the two-mass-easy example, we show that overall momentum conservation is 
an important constraint on the solutions of spinor components. In the two-mass-hard ex- 
ample, we show that any pair of adjacent three-gluon corners must come with opposite 

helicity configurations (one each of + H — and h), which is consistent with Figure 4. 

This restriction figures into our one- mass example also, where we choose to compute a box 
integral whose quadruple cut has a next-to-MHV amplitude at one corner. 
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The four-mass coefficient presented here is new, but our examples in this section for 
three-, two-, and one-mass coefficients are all for a next-to-MHV seven-gluon amplitude that 
was computed in f|] and f|. Our formulas reproduce the ones in those papers exactly. 

In section 4, we will present several examples of infinite classes of coefficients for 
n-gluon next-to-MHV amplitudes. 

3.1. Four-Mass Box Integral Coefficients 

Let us illustrate our solution by computing the coefficients fi , / 2 of the scalar box 
functions /f^^-i an d ^2:2:2-2 f° r t ne helicity configuration (1~, 2~, 3~, 4~, 5 + , 6 + , 7 + , 8 + ). 
We can see immediately that f± = 0, because any possible assignment of helicities for the 
internal lines, one of the four tree-level amplitudes in ( |3.1| ) always vanishes. 

The box diagram for Ig™-2-2-2 * s shown in Figure 5. 




Fig. 5: The box function /f™ :2:2 ;2 for the helicity configuration (1", 2~ , 3 , 4", 5 + , 6 + , 7 + , 8+) 
with internal momenta ^1,^2,^3,^4, and its associated quadruple cut. 



This helicity configuration allows only gluons to circulate in the loop. Moreover, each 
vertex is MHV, so it is straightforward to write the coefficient ( |3.1| ) in terms of Parke- Taylor 
formulas [ED . 



/ 2 = ^ tree ((-£ 4 )+,8+,i-,^)x 

x A tree ((-£ 1 ) + ,2-,3-,(-£ 2 )+)A tree (£ 2 -,4-,5+,(-£ 3 ) + )A tree (£3,6+,7+,£ 4 - 
(Hi) 3 (2 3> 3 (£ 2 4> 3 (£ 4 £ 3 ) 3 



= -y - 

2 (h 8) (8 1) (3 £ 2 )(£ 2 2) (4 5) (5 £ 3 )(£ 3 £ 2 ) (£ 3 6) (6 7) (7 £ 4 ) 

= 1 V [6 7] 3 (1|A / 2 |4) 3 

2 2^ (8 1)[2 3] (4 5)<5|/ 3 k k\W>\h / 3 |6] <8|/ 4 |7] ' 

(3.4) 
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Further information about four-mass box integrals, consistency checks for this coeffi- 
cient, and twistor space structure may be found in the appendices. 

3.2. Three-Mass Example 

Consider the three-mass box scalar integral for the seven-gluon next-to-MHV ampli- 
tude A 7: i(l-,2-,3-,4 + ,5 + ,6 + ,7+) shown in Figure 6. 




Fig. 6: A three-mass scalar box integral for the next-to-MHV amplitude 

A 7:1 (l-,2- 3-,4+,5+ 6+ 7+). 



As explained in section 2, for the three-gluon amplitude with helicities (+ H — ),we 
should choose the A's to be proportional so that the amplitude does not vanish. We make 
this choice in solving the equations p 2 ■ £\ = p 2 • £4 = £\ • £4 = 0, which can be written as 
(h 2}[£ t 2] = (£ 4 2}[£ 4 2] = (£ t £ A )[£ X £ 4 ] = 0. 

We choose to write each four gluon amplitude in the generalized unitarity cut diagram 
as a mostly minus MHV amplitude. This gives 

~ 1 ( [*i 4] 3 \f [4 4] 3 \ ( [5 6] 3 \( [£3 7f 



[h 2] [2 £ 4 ] J V [I2 h] [h 3] [3 4] ; V [6 4] [4 h\ [£2 5] J\ [7 1] [1 £4] [£4 £3] 

(3-5) 

Now we explain how to solve \i 1 and A^ from light-cone conditions. First, using 
= 2£ x ■ p 2 = (2 £ x )[2 h], we get (2 £ x ) = or [2 £ x ] = 0. To get a nonzero contribution, 
we choose to solve (2 £1) = 0, so that = a\ 2 . Now we can solve for a: 

2£x ■ (P3+P4) = (P3+P4) 2 = -a(2|(3 + 4)|4]; 

-(P3+P4) 2 (3-6) 



a = 



(21(3 + 4)^] 
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To solve for A^ , we use the equation £\ = (£i — (ps + P4 + p$ + Pe)) 2 = 0, i.e., 

2h ■ (p 3 +P4+P5+ Pe) = (Pa +P4+P5+ Pg) 2 = -a(2\(3 + 4 + 5 + 6) \h] 
(ps +P4+P5+ Pe) 2 (2|(3 + 4) |4] = (p 3 + p 4 ) 2 (2| (3 + 4 + 5 + 6) \£ x ] 

(3.7) 

Note that this is not enough to fix the whole Xi 1 , but it fixes its direction. This is good 
enough to compute (|3.5|), since the coefficient has degree zero in A^. 

Then we do the same for £4 and compute A^ 4 and \i 4 up to a scale. 

After this is done, we choose to write 

_ (2|l 2 |.] _ (2J4M , , 

[£2 ] - 724P [4 #] - 724T- (3 - 8) 

where • represents any external gluon, £\ or £4. It is easy to see that all denominators will 
drop out from the fact that the coefficient has degree zero in £2 and £3. 
This gives for the coefficient ( |3.5|) the following: 

~ _ 1 ( ['1 4] 3 \ I <2|^2|4] 3 \f [5 6] 3 \ / (2|£ 3 |7] 3 

93 



2 V[4 2] [2 4]; \ v (2|£ 2 |£ 1 ][£ 1 3][3 4]y V (2|4|6](2|4 ■ 4|2)(2|£ 2 |5] ; \J7 1][1 4](2|£ 3 |4] 

(3.9) 

where we will substitute for £2 and £3 using 

4=4-^3-^4, *3=4+Pl+P7- (3-10) 

The coefficient of the associated three-mass scalar box function was computed by two 
different methods in [|]||, and it is given by 

£3 

93 = 



(P3+P4+P5 +Pg) 2 (P2+P3 +P4) 2 - (P3+P4) 2 (P7 +Pl) 2 

(1 2) 3 (2 3) 3 [5 6] 3 

(7 1)(3 4)(2|3 + 4|5](2|7 + 1|6]((7 1) (2|3 + 4|1] - tf (7 2))(t 3 4] (2 4) - (3 4)(2|7+ 1 1 3] ) 

(3.11) 

It is easy to check that the two formulas ( |3.9|) ,( |3TTl~l) give exactly the same answer for 



93- 



3.3. Two-Mass-Easy Example 

Consider the coefficient C4 in the amplitude ^7.1(1", 2~, 3~, 4 + , 5 + , 6 + , 7 + ). The box 
integral is shown in Figure 7. 
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(a) (b) 

Fig. 7: Quadruple cut of a two-mass-easy scalar box integral for the amplitude 

A 7:1 (l-,2-,3-,4 + ,5+,6 + ,7+). 



Notice that there are two possibly nonvanishing helicity assignments for the cut prop- 
agators. Here it is possible to solve for the spinor components of (^1,^2,^3,^4), up to a 
scale for each, for both helicity configurations. However, the solution for configuration 
(b) found from this procedure fails to satisfy momentum conservation at the "massive" 
corners. Therefore, there is only one true solution, namely configuration (a). From it we 
find 



~ _ 1 ( (h h) 3 \ ( [h 6] 3 ^ ( [i 3 7] 3 ^ ( [h £ 4 

C 4 — 77 



3 



2 V (h 4) (4 5) (5 £ 2 ) J V [6 £3] [h h] ) \[7 1] [1 2] [2 £ A ] [l A £ 3 ] ) \ [£ 4 3] [3 £ x \ 
Then 



(3.12) 



' (P4 +P5 +Pe) 2 {P3 +PA +P5) 2 - (P4 +P5) 2 {P7 + Pl +P2) 2 ' 

which agrees with the results given in B,B| . 

3.4- Two-Mass-Hard Example 

Consider the coefficient d^^ in the amplitude A 7: i(l _ , 2~, 3~, 4 + , 5 + , 6 + , 7 + ). The 
box integral is shown in Figure 8. 

There are two possibly nonvanishing helicity assignments for the cut propagators. For 
configuration (b), it is not possible to solve for the spinor components of (^1,^2,^3,^4)- 
The reason is that from the pi corner, we derive the condition that \i 3 , Xi 4 and A7 are 
all proportional, and similarly from the p\ corner, we derive that A^ 4 ,A^ and Ai are all 
proportional. But this means that A7 is proportional to Ai, which is not true of generic 
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(a) (b) 

Fig. 8: Quadruple cut of a two-mass-hard scalar box integral for the amplitude 

A 7:1 (l-,2-,3-,4+,5 + ,6 + ,7+). 



external momenta. The lesson here is that adjacent three-gluon box corners must have 
opposite helicity types in order to give a nonvanishing product of amplitudes. 

For the helicity configuration (a), there is no such obstacle to solving for the spinor 
components of (^1,^2^37^4) up to a scale, because the conditions are that \i 3 ,\i 4 and 
A7 are all proportional, but it is the antiholomorphic spinors A^ 4 ,A^ 1 and Ai that are 
proportional for the other corner. Once the solution is obtained, the coefficient ds y2 is 
given by 

1 ( (h I) 3 \ ( \h 7] 3 \ / (£3 £2? \f (2 3) 3 



d 



3 ' 2 2 \{lh)(hh)J \[hU}[hii]J V(^5}(5 6)(6 4)7 \(3 4)(4£ 2 )(£ 2 £ 1 )(£ 1 2) 

(3.14) 

Then 

da,2 = —2-, ; ; — 3 -%t ; ttt, (3.15) 

\Ps +P6 +P7) (P7 + P1) 

which agrees with the results of 

3.5. One-Mass Example 

Even for the one-mass box functions, where the internal momenta might seem to be 
over constrained, we can solve for the coefficients using generalized unitarity. Consider the 
coefficient 62 in the amplitude At-a(1~ , 2~, 3~, 4 + , 5 + , 6 + , 7 + ). The box integral is shown 
in Figure 9. We have applied the lesson learned in the previous subsection, that adjacent 
three-gluon box corners must have opposite helicity types in order to give a nonvanishing 
product of amplitudes. Thus only the two helicity configurations shown in the diagram 
might contribute to the coefficient. 
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(a) (b) 

Fig. 9: Quadruple cut of a one-mass scalar box integral for the amplitude 

A 7:1 (l-,2-,3-,4+,5 + ,6 + ,7+). 



For the configuration (a), we solve for the spinor components of (^1,^2,^3,^4) while 
requiring that 

[h 6] = [is 6] = [£ 2 £3} = 0, (£3 7) = (£ 4 7) = (£ 3 £4), [£4 1] = [h 1] = [£4 £1] = 0. 

(3.16) 

Then we substitute this solution into the product of tree amplitudes: 

«-) _ 1 ( ( 2 3 ) 3 ^ f tela) 3 N ( ih7? \ ( (Hi) 3 



J ' 

2 2 \ (3 4)(4 5)(5£ 2 )(£ 2 £ 1 )(£ 1 2)J \(£ 2 6)(6£ 3 )J \ [7 £4} [£4 £3} J \(£iii)(hl) 

(3.17) 

For the configuration (b), we solve for the spinor components of (^1,^2,^3,^4) while 
requiring that 

(£2 6) = (£3 6) = (£ 2 £3) = 0, [£3 7] = [£4 7] = [£3 £4}, (£4 1) = (£1 1) = (£4 £1) = 0. 

(3.18) 

We use this solution in the product of tree amplitudes, which now includes an NMHV 
amplitude [p2| , |23| : 

/3 2 



4 ree (Mi)",2-,3-,4+,5+,£+) 



t 5£ 2 (-l 1 ) S U 2 s e.2{-tx ) s 23«34 



te 2 (-£ 1 )2Se 2 (-£ 1 )S(-e 1 ) 2 S34S4 5 S45S 5 £ 2 S£ 2 (_^ 1 )S(_^ 1 ) 2 S23S34 

0=[5 £ 2 }(2 3)(£ 1 \5 + £ 2 \4], 
7= [4 h){£ x 2>(3|4 + 5|4], 

Ujk = (i 3} + (i k)[i k] + (j k)[j k}. 

(3.19) 
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The contribution of this configuration to the coefficient is 



-1 



(3.20) 



6 2 a) +6 2 b) 



(P6 +P7) 2 (P7 +Pl) 2 , 



Then 



which agrees with the results of We note that in [|J, this coefficient 62 was given in 

terms of a seven-gluon tree amplitude and eight other coefficients (using a relation derived 
from the infrared singular behavior of a certain unitarity cut). Here, the quadruple cut 
allows us to obtain an explicit formula for 62 directly. 

One can easily obtain 62 for the n-gluon next-to-MHV amplitude (1~, 2~, 3~, 4 + , n + ) 
by substituting the tree amplitude A%l\((-£ 1 )-,2-,3-,4+, (n - 2)+,£+) for A^ ee in 
( p.20|) . Compact expressions for that amplitude have been given in [|l4"l ,|il 



4. All-Multiplicity Examples for NMHV Amplitudes 

In this section we present some examples of using quadruple cuts to compute some 
classes of coefficients for n-gluon next-to-MHV amplitudes. Here we substitute the actual 
solutions for the cut propagators, so that our final formulas are given in terms of the 
external momenta only. 

4-1. All- Multiplicity Examples of Three-Mass Coefficients 

In this subsection we present three classes of three-mass coefficients for next-to-MHV 
n-gluon amplitudes. The first two are depicted in Figure 10. 




Fig. 10: Two families of three-mass scalar box integrals in an NMHV configura- 
tion. 
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Consider first the three-mass scalar box integral coefficients g n :r:r';i-, where the three 
negative helicities are j, k, i — 1, and j is in the vertex of (i, + r — 1) and k is in the 
vertex of (i + r + r', i — 2). See Figure 10(a). Using (|3.1| ), we can write down 

(i - 2 % - - 1 + r - 1 % + r)(i + r + r' - 1 i + r + r') [£i £4 3 
(h J') 4 (£3 £2? (k 



fn:r:r':i 



X 



\£2 £i)(£i i)(i + r-l£ 2 ) (£2 i + r)(i + r + r' -1 £ 3 ) (£ 3 i + r + r')(i - 2 h)(£ A £3)' 

(4.1) 

We simplify the above formula by noticing that 

(p\K 2 ■K 1 \i-1) = (p £ 3 ){£ 3 £ 2 ](£ 2 1 - 1) 

(4.2) 

(p\K 2 -K 3 \i-l) = (p £ 2 ){£ 2 £ 3 }(£ 3 i-l) 
and X £l = aAj_i, \t 4 = (3\ i - 1 , so 

(i - 2 i - - 1 + r - 1 % + r)(i + r + r' - 1 % + r + r') (i - 1 j) 4 (k i - l) 4 



in:r:r' \i 



m n s =x(ss + l) (»-li>(»-2t-l> 

(-KI) 3 



(i + r - 1|K 2 ■ i^sN -!><* + r \K 2 ■ K 3 \i - 1) 

1 

X (i + r + r'-l\K 2 -K 1 \i-l)(i + r + r'\K 2 ■ K x \i - 1) 

x ^i¥ L Ty (aV(I - 1|4 '^- 1>) ' 

(4.3) 

We will show below that the last line is simply {K\ X K\ 2 — KfK^j. Thus we get immedi- 
ately the coefficients 

j) 4 (k i - l) 4 (z + r-li+r)(i + r + r'-li + r + r')(K^) 3 
= n:=i(^ + l>(i + r-l|^2-^3|i-l> (44) 

X (i + r\K 2 ■ K 3 \i - + r + ? - 1\K 2 ■ K x \i - + r + r'\K 2 ■ K x \i - 1) ' 

It can be shown that if we put j = 3, i — 1 = 2, k = 1 we reproduce the all-multiplicity 
NMHV coefficients given in ||. Also, our general formula matches the results of for 
seven gluons, i.e., the coefficient of helicity assignment (1~, 2~, 3 + , 4~, 5 + , 6 + , 7 + ) and 
C257 of helicity assignment (1~, 2 + , 3~, 4 + , 5~, 6 + , 7 + ). 

Now we prove the identity involving the last line of ( |4.4| ), namely that 

[h ■ i £i] (« 2 /3 2 <* - 1|4 • < 2 |* - 1)) = KlK* 2 - K 1 2 X 3 2 - (4-5) 
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To do this, we need to notice that 



1~ a~ 

\e 4 = p\i-i+p\ ei , (4.6) 

which can be obtained by noticing that £4 = £\ + Pi-\ while Pi-\ = Ai_iAi_i, £\ = 
aXi-iX^, £4 = f3\i-i\e 4 . Using this we can easily show that 

(K 2 41 K^ - K\K\) = K\(2 Vl _ x ■ K 3 ) + Iff (2 Pi _x ■ K x ) + {2 Pi _ x ■ K z ){2 Pi _ x ■ K x ) 

= (-) {h h) [£1 £2} {i - 1 4) [i ~ 1 h] + (-) (h £3) [U £ 3 ] (-) {i - 1 h) [i - 1 £2] 
+ (i-l£ 3 )[i-l£ 3 ](-)(i-l£ 2 )[i-l£ 2 ] 

= -a(i-l £ 2 )[£i £2] (i - 1 4) [i ~ 1 4] + 0<i - 1 4) [4 4] (i - 1 £2) [i - 1 £2] 

- £ 3 )[i-l £ 3 }{i-l £ 2 )[i-l £2] 

= -a(i -l£ 2 ){i-l £ 3 )[h i - 1][£ 2 £ 3 ] =a[i-l - 1 £ 3 )[£ 3 £ 2 }{i - 1 £ 2 ) 



(4.7) 



and 

[ilh]3 (a 2 (3 2 (t-l\£ 3 -£ 2 \t-l)) 



[£ 4 t-l][i-l£i] 

= a[z-l £i](z-l £ 3 )[4 

By similar calculations, we have found another series of coefficients of F 3m . where the 
three negative helicities j, k,i — l are at the vertices (i+r, i+r+r' — 1), (z + r+r', — 2) 
and massless leg (i — 1) respectively. The coefficient is given by 

_(k i - l) 4 (j\K 2 ■ Ki\i - l) 4 (i + r - 1 i + r) (i + r + r' - I i + r + r') 
9 ™ r "> i ~ (-K!)U n s=1 (ss + l) 

x \ (4 9) 

(i + r -1\K 2 - K 3 \i-l)(i + r\K 2 - K 3 \i-1) V ' ' 

\ 

{i + r + r' - l\K 2 ■ K\\i - + r + r'\K 2 ■ K x \i - I)' 

Now we offer an example of box integrals whose quadruple cut involves fermions and 
scalars, in addition to gluons, circulating in the loop. 
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(a) 



(b) 



Fig. 11: Another family of three-mass scalar box integrals in an NMHV configu- 
ration. 



In this case, the three negative helicities j, k,n are at vertex K\,K 2 and K 3 respec- 
tively. There are two possible helicity configurations for loop momentum. See Figure 11. 
Using Ward identities, one can prove the following relations between amplitudes involving 
gluons and amplitudes involving fermions and scalar s p4,^3|. 



A(Fr,g+,..,gT,...,F+) 
A[Si , g 2 , ••) gj i •••) S£) 



(j n) 
O'l) 
(j n) 2 
(J I) 2 



A M ^(g^,gt^..,gr,... : gt), 
A M ^(g^g+,...,g-,...,g+) 



(4.10) 



we can write down the sum of configurations (a) and (b) as 
(i — 2 i — 1) (i — 1 i) (i + r — 1 i + r) (i + r + r' — 1 i + r + r') 

2 nr=i (**+i> 

(hj) 2 (£ 2 j) 2 (£ 2 k) 2 (£ 3 k) 2 {£ 3 n) 2 (£ 4 n) 2 



x 



i)(i + r-l £ 2 ) {£2 £1) {£2 i + r) (i + r + r' - 1 £ 3 ) (£ 3 £ 2 ) (£ 3 i + r + r') (i-2 £ 4 ) (£ 4 £3) 



[£ 4 i - l] 2 [h i - l] 2 
£!}[£! £4] [U 

where 



xA~ 2 (A-l)\ 



(4.11) 



A 



(j^jk £ 2 )(n £ 3 )[i-l£i 
(j £ 2 )(k £ 3 )(n £ 4 )[i-l£ 4 



(k\K 2 -K 3 \i 


-l){n\K 2 -K x \i 


-1) 


(j\K 2 -K 3 \i 


-l){k\K 2 -K x \i 


-1) 
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(4.12) 



After simplifying the above expression similarly to the previous examples, we obtain 

Qn:r:r':i 



({k\K 2 ■ K 3 




i - i 


-l)-{j\K 2 -K 3 


i-l){k\K 2 -K 1 


% - 


-l)(nz-l)) 4 


(i + r-l\K 2 -K 3 


i ■ 


-l)(i + r 


K 2 -K 3 


i-l)(i + r + r' -1\K 2 -K X 


i-1) 



X 



(i + r — 1 i + r) (i + r + r' — 1 i + r + r') 



(i + r + r'\K 2 -K^i- l)K*(i - 1\K 2 ■ K 3 \i - l) 4 n"=i<^ s + 1) 



(4.13) 

To simplify further, we use the following result: 
(k\K 2 ■ K 3 \i - l)(n\K 2 -K^i- (j\K 2 ■ K 3 \i - l)(k\K 2 ■ K x \i - l)(n i - 1) 

= -(i - 1\K 2 ■ K 3 \i - 1) ((k\K 2 ■ K 3 \n)(j i - 1) + (k\K 2 ■ K x \j){n i - 1) 

-(ni-l)(j i-l)(k\K 2 \i-l}). 



(4.14) 



We arrive at the final expression 

Qn:r:r':i 



({k\K 2 -K 3 


\n)(j + (k\K 2 ■ (Kx +p i - 1 )\j){n i 


- 1)) 4 




(i + r - 1\K 2 ■ K 3 \ 


i — + r 


K 2 -K 3 


i - + r + r' - 1\K 2 ■ K\ 


i-i) 



X 



(i + r — 1 i + r) (i + r + r' — 1 i + r + r') 
(i + r + r'\K 2 ■ Kt\i - l)KlYT s=1 (s s + 1) 



(4.15) 

We have checked that this expression reproduces the coefficient C135 of seven gluons in the 
helicity configuration (1~, 2~, 3 + ,4~, 5 + , 6 + , 7 + ) as given in |J. 

4-2. An All- Multiplicity Example of a Two-Mass-Hard Coefficient 

Now we consider a series of F 2m h shown in Figure 12, where both massless legs (i — 1) 
and (i — 2) are negative helicities. 




i-2 K; 



K- 



i-2 K; 



+ 




+ 



+ 




+ 




i-1 



K 4 



(a) 



(b) 



Fig. 12: Quadruple cuts of a family of two-mass-hard scalar box integrals in an 
NMHV configuration. 
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The third negative helicity j is at the vertex of (i, % + r — 1). For this case there 
are two contributions. For part (a) we have 

_(i-3i- 2)(i -2i- l)(i - 1 i)(i + r- li + r) 

\hUf (hj)* (hh) 3 (i-2£^ 



x 



[£ 4 i - 1] [i - 1 h] (h i)(i + r-l £ 2 ) (£ 2 h) (h t + r)(i-3 £ 3 ) (£ 3 i - 2) (£ 4 £ 3 ) ' 

(4.16) 

Using similar manipulations, with = cuAj-i, Xi 4 = (3\i-x and the identity Q15|) , which 
here takes the form 



Vl i-W-l tl] X ^ ~ 1 ^ ^ ' ~ 1) = A « A '« " KlKl (4.17) 

(K-l+iflf^+Xj) 2 



we can read the coefficient from part (a) as 

(K 2 ) 3 (i-3 i — 2> <i — 2 z-l) 4 (z-l j) 4 (i + r-li + r)/lT s=l (s s + 1) 

(i + r- 1\K 2 ■ Pi _ 2 \i - l)(i + r\K 2 ■ p t _ 2 \i 3\K 2 ■ K x \i - l)(i - 2\K 2 ■ K x \i - 1) 

(4.18) 

For part (b) we must use different spinor relations, Xi 4 = 5A;_2, Xi 3 = fi\- 2 and an 
analog of the identity (|4.5|) , namely 



[<4 ' :jJ " :a 2 (3 2 (t - 2 £ 2 )[£ 2 i-2) = K 2 ,K 2 X - KlKl 



[f 3 i-2][i-24] rr -^^x-x. -/ "4^, (4ig) 

= (Pi-2 +Pi-i) 2 (Pi-i + Kx) 2 , 

to read out the coefficient 





|z-2) 4 (z-2 i- 


-l) 4 (i-l i)(i + r-l i + r)/(K 2 YT s=1 (s 8 + 1)) 


(i-l\K x -K 2 \ 


1-2)^! -K 2 \ 


i-2)(i + r-l\K 1 -p i - 1 \ 


i-2)(i + r\K 1 -p l _ 1 \ 


i - 2) 



(4.20) 

Putting it all together, we finally have 



dn:r:i 



(Kj) s {i - 3 i - 2)(i - 2 i - l) 4 (i - 1 j) 4 (i + r - 1 i + r) / Us=i( s s + x ) 
(i + r- 1\K 2 - Pl _ 2 \i - l){i + r\K 2 -pi-2\i -l)(i- 3\K 2 ■ K x \i - l)(i - 2\K 2 ■ K x \i - 1) 

+ 



(j\KfK 2 


\i - 2) 4 (z -2 %- 


- l) 4 (i -li)(i + r-li + r)/(K 2 Us=i( s s + 1)) 


(i-l\K x -K 2 \ 


i-2)(i\K 1 -K 2 \ 


i-2)(i + r-l\K 1 -p i _ 1 \ 


i-2)(i + r\K 1 -p i _ 1 \ 


i - 2) 



(4.21) 

We have compared this general formula with coefficient C457 of seven gluons in the helicity 
configuration (l _ 2 _ 3 + 4~5 + 6 + 7 + ) in and found that they agree. 
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5. Summary 

In this paper we have reduced the problem of computing the coefficient of any scalar 
box integral in any one-loop Af = 4 amplitude to finding solutions to the four equations 

£ 2 = 0, (£-ifi) 2 = 0, (£ - K x - K 2 f = 0, (£ + K 4 ) 2 = (5.1) 

in ( h+) signature. From this set of solutions it is possible to read off the coefficient 

from the formula 

^ = j4[ E njAT e AT e A^ e AT e . (5.2) 
' ' S,J 

It would be interesting to perform a full classification of helicity configurations to obtain 
explicit formulas for all coefficients to all multiplicities. 
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Appendix A. Discontinuities of the Four-Mass Scalar Box Integral 

In this section we study in detail the discontinuities of the four-mass box integral. 
The motivation is to allow us to perform non-trivial consistency checks on the coefficients 
found in the previous section and to discuss the difficulties in trying to compute them by 
the application of coplanar operators along the lines of |||4[]. 

One way to obtain the discontinuities is to compute the imaginary part of the explicit 
formulas for the four-mass integral in a kinematical regime chosen in order to isolate the 
cut of interest. Another way, which we find more intuitive, is to compute the discontinuity 
of the scalar box integral by cutting it directly. It turns out that the cut integral is easy 
to evaluate explicitly, as we now describe. 
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A.l. Cut Integral 

Consider first the four-mass box function, 



r4m 



dH (£ 2 + ie){{£ - K x y + ie)((£ -K x - K 2 ) 2 + ie)((£ + K 4 ) 2 + ie) " (A-1) 



Let us start with the discontinuity in the (K\, i^-channel. The integral we have to 
evaluate is obtained from (|A.1| ) by replacing the first and third propagators by the delta 
functions imposing the on-shell condition. 



\/-><"l / V (+ W + H(l-ifl-if2) 2 ) 



r4m I / j4 



Now we can parameterize £ a a = tX a Xa and write the Lorentz invariant measure as 
follows: 

r°° f ~ 5^((£ - K-i - Ko) 2 ) 

n?,>° = i t dt J (\ d\)[\ dX] (< \ i)2( ; + 4V - (A.3) 

The remaining delta function can be written as 

5((£ -K x - K 2 f) = 8{tX a X h K™ - K 2 12 ) = -J 7 8{t - 5 12 . ), (A.4) 

A a A a i\ 12 A a A a J\ 12 

where Ki 2 = K\ + K 2 . Therefore the t integral can be performed to find 

A/ 4m k2 >0 = f(X dX)[X dX] — = ^ . (A.5) 

12 J \X a X & K™) 2 {£-K x ) 2 {£ + K A ) 2 

The denominator can be written as 

(X a X d K™) 2 (£ - i^) 2 (£ + K 4 f = (Q ad X a X d )(S aa X a X a ), (A.6) 

with 



Qaa = _ K f 2 Kf + K 2 K™, 

S ad = K\ 2 Kf + K\Kl%. 
The way to evaluate this integral is to use a Feynman parametrization. This gives 



(A.7) 



A/ 4m | K 2 >0 = f da 1 da 2 6(a 1 + a 2 - 1) I (X dX) [X dX] ^ — , (A.8) 

Jo J (P ad X a X d ) 2 

with 

P ah = ai Q ad + a 2 S ad . (A.9) 
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The integral over the sphere was evaluated in [ 13 1 and gives 1/P 2 . Therefore we have 



f 1 1 

A/4m |K? 2 >o= / <faitfa2<y(tti + Q2-l) 2n2^nn o T 2 Q2 - ( A - 10 ) 

Now we can perform the integration in a± and a 2 . The result is 

r4m, 1 -A 2 -p) 2 



where p is defined as in Ql.4|) . We can similarly evaluate the discontinuity of the four- mass 
box function in the Xi-channel. This is obtained from ( |A.1|) by replacing the first and 
second propagators by the delta functions imposing the on-shell condition: 

AI \k>« = j dl { i- Kl -K 2W + K,y (A - 12) 

Again, parametrize £ aa = tX a X a and then perform the t integral to find 

AI 4m \ K 2 >0 = / (A dX) [A dX] — = x - . (A.13) 

Here the denominator can be written as 

(XaKKffii-Kx -K 2 ) 2 (£ + K 4 ) 2 = (Q aA X a Xa)(S a& X a Xa), (A.14) 

with 

Q ah = K 2 2 Kf - K 2 K™, 

-2 jsaa i 7^2 jraa 



(A.15) 



gaa = K ^ K aa + K Z R a, 

So we may follow exactly the same steps as in the previous case, and here we find the 
result 

A r4m, , _ 1 Inr l^ (1 ~ Al + A2 ~ p)2 ^ ( \ 

AI lK >°--^KlKT 2 l ° g { 4X 2 J' (A " 16) 

where Ai, A 2 ,p are the same expressions defined in (P-.4Q . 

A.2. Imaginary Part 



As a consistency check, we can reproduce our results (|A.11| ) and ([A.16Q from the 



imaginary part of the explicit form of the four-mass integral in the appropriate kinematical 
regimes. 
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The integral (|A.1| ) was computed explicitly in [25] and is given in terms of the invari- 
ants t m i = — {K m + K m+1 + . . . + Ki-i) 2 as follows: 



r4m 



1 



E(-l) j {-\\n 2 [-x s ) 



Li 2 ( H -x 

ti3 - ze 

Ll 2 1 + ~Xj 

ti2 - ie 



Ha - ie \ ( t 34 - ie 
j I V \ ~ x k, ~. I In ( 1 + -x, 



V -x k , 



ti3 - ^ 
t 2 4 - ie 
h2 - ie 



ti3 ~ ie 



In ( 1 H — — -x- 

ti2 - ie J 



ax 2 + bx + c + ied = a(x — x\)(x — X2), 



with 



(A.17) 



+ ln(-Xj-)(ln(i 1 2 - ie) + ln(t 13 - ie) - ln(t M - ie) - ln(t 23 - ie))) . 
Here the function rj(x,y) is given by 

r)(x,y) = 2jvi['d(—Imx)'d(—Imy)d(Im(xy)) — i?(Ima;)i?(Imy)i?(- lm(xy))], (A. 18) 

and x\ and x 2 are the roots of a quadratic polynomial: 



(A.19) 



(A.20) 



a = £24^34, 

b = ^13^24 + ^12^34 — ^14^23) 
C = ti2^13j 
d = t 23 . 

The ie prescription in (|A.17| ) allows us to use this formula in any kinematical regime. 



The main simplification is that the proper branch of each of the functions in (|A.17|) is simply 
given by the principal branch. The formula for I Am presented in || can be recovered from 
( |A.17p by setting the r\ functions and e to zero. Even though the formula looks simpler 
in that form, it has the disadvantage that it is not clear which branch of the various 
dilogarithms and logarithms compute the appropriate discontinuity. 



Appendix B. Twistor Space Structure of Four-Mass Box Coefficients 

In this section we study the twistor space localization of the four-mass scalar box 
function coefficient found in section 3.1. The way we choose to do this can also be thought 
of as a consistency check on the coefficient. 
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Consider for example the amplitude A 8; i(l~, 2~, 3~, 4~, 5 + , 6 + , 7 + , 8 + ). The two four- 
mass box integrals have momenta distributed as (12) — (34) — (56) — (78) and (23) — (45) — 
(67) — (81). In the first case we find that the coefficient must be zero, since one of the 
tree-level amplitudes in the quadruple cut diagram is necessarily zero. This is essentially 
equivalent to the observation of M, where the same conclusion was derived from a triple 
cut. 

In order to study the second four-mass integral, consider the cut in the (2345)-channel. 
This cut is given by 

C 23 4 5 = J dfi A t 6 ree (£ 1 ,2-,3-,4-,5+,£ 2 )4 ree (-£ 2 ,6+,7+,8+,l-,-£ 1 ). (B.l) 

There are four contributions to this cut depending on the helicities (h^, hg 2 ) of the 
particles running in the cut propagators. If {h^.h^) = (— , — ) then we get zero. If 
(h^jhtz) = (— , +) or (h^jhiz) = (+,—), then the whole Af = 4 supermultiplet con- 
tributes. In this case, both tree-level amplitudes are very simple; they become a mostly- 
minus MHV and a mostly-plus MHV, respectively. These contributions can be easily 
written as 

(l|(2 + 3 + 4)|5] 4 

Since the four- mass integral does not contribute to the cut C1234, then it does not 
contribute to this part of the cut C2345 either. 

Finally, we have the case with {h^.h^) = (+,+). In this case, only gluons can 
propagate; the complication arises from the fact that both tree-level amplitudes are next- 
to MHV six-gluon amplitudes. 

It turns out that we cannot use collinear operators to extract information from this 
cut. The reason is that a collinear operator does not localize the integral. We need a 
coplanar operator. 

A coplanar operator is defined as follows O]. 



K ijkl = (i j)[d k di] + (j k)[di d^ + ik i)[dj di} + (k l)[di djj + ii l)[dj d k ]+(j l)[d k 8,1 (B.3) 
where 

(*)* = ^ 

In this case we could use ^2345 or i^678i to produce rational functions. 
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Recall that the cut can also be written as the discontinuity of the amplitude in the 
(2345) -channel. 

C 23 45 = AAl- loop = ... + fAI 4m . (B.5) 

We have only explicitly written the term from the four-mass integral. Once we apply either 
of the two coplanar operators K on the cut integral, we produce a rational function. On 
the other hand, acting with K on (|B.5|) produces 

KC 2345 = ... + K(fAI 4m ). (B.6) 

Following the arguments of , one can show that the terms in the ellipses, which come 
from lm, 2m and 3m scalar box integrals, are rational functions. Therefore we conclude 
that 

K(fAI 4m ) (B.7) 

must be rational. 

But we can go even farther. In the previous section, we found that 



A r4m _ 1 i ( (1 ~ Aj ~ A 2 ~ pY , 

A/ ~W!M, g \ ^ '■ ( 

Therefore, when none of the derivatives in K act on the logarithm, we find a term of the 
form 



/ 1 wl /(l-Ai-A 2 -p) 



P K{ X K{ 2 ) ° V 4AiA 



K .J y2 x log ^ L - / " . (B.9) 



\2 



The only way this is consistent with ( |B.7| ) being rational is that 



f_ 



K ( -^TITT 1 = 0- (B.10) 



Recall that / = —2f/(pK2 1 Kf 2 ) is the definition of the four-mass box function co- 
efficient, which we claimed would have a simple twistor structure configuration. Indeed, 
( [B.10Q is the reason for our claim. 

So we have found that -^2345/2 = -^6781/2 = 0. Note that this four-mass box integral 
also has a cut in the (4567)— channel. Using the same logic we find that -^4567/2 = 
-^8123/2 = 0. 

The conclusion is then that all gluons in two adjacent corners of the four-mass box 
function coefficient are localized on a plane. The most general configuration consistent 
with this picture is shown in Figure 13. 
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Fig. 13: The twistor space structure of the four-mass coefficient fi- The four lines 
are not all coplanar. 

The consistency checks we have run on the coefficient found in section 3.1 are the 
following: f 2 is annihilated by the coplanar operators mentioned before, and ( |B.7| ) is a 
rational function. 

At this point we must comment on why studying this cut does not give an efficient 
way of computing f 2 . After all, computing the rational function -^2345^345 from the cut 
integral representation of C2345 is very simple. Moreover, it turns out that all other coeffi- 
cients contributing to this cut can be computed from cuts in three-particle channels, along 
with linear equations from the infrared singularities, using the method of ||. Therefore, 
all reduces to 

W = K(fAI 4m ) : (B.ll) 

where W is a known rational function, and AI 4m is also known. 

Let us write K schematically as 0\0 2 - Then, expanding the right hand side of ( [B.ll] ) 
we find 

W = K(f)AI 4m + 1 (/)0 2 (AI 4m ) + 2 (f)0 1 {AI 4m ) + fK{AI 4m ). (B.12) 

The problem arises when one notices that the only term that has no derivatives acting 
on / vanishes. Indeed, it not difficult to show with the help of a symbolic manipulation 
program that K(AI 4m ) = 0. 
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